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The small t asymptotics of a class of solutions to the 2D cylindrical Toda 
equations is computed. The solutions, qk(t), have the representation 



where are integral operators. This class includes the re-periodic cylindrical 
Toda equations. For re = 2 our results reduce to the previously computed 
asymptotics of the 2D radial sinh-Gordon equation and for n = 3 (and with 
an additional symmetry constraint) they reduce to earlier results for the radial 
Bullough-Dodd equation. Both of these special cases are examples of Painleve 
III and have arisen in various applications. The asymptotics of qk(t) are derived 
by computing the small t asymptotics 



where explicit formulas are given for the quantities a k and b^. The method 
consists of showing that the resolvent operator of K k has an approximation 
in terms of resolvents of certain Wiener-Hopf operators, for which there are 
explicit integral formulas. 

Introduction 

We consider here solutions of the cylindrical Toda equations 



Abstract 



q k (t) = log det (I - A K k ) - log det (I - A K k ^) 




ql(t)+t- 1 q> k (t)=A(e 




(i. 



satisfying the periodicity conditions qk+ n = Qk- The integer n is arbitrary but fixed. 
It follows from results in || that solutions valid for all t > are given by 

q k (t) = log det (/ -XK k )- log det (/ - \K k - X ), (1-2) 
where K k is the integral operator on R + with kernel 

$>*^ £ , (1.3) 

u -uju + v 

u running over the nth roots of unity other than 1 . 

In the case n = 2 we have q k +i = —qk and ( |1.1| ) becomes, with q equal to either 

Qk, 

q"{t) +r 1 q'(t) = 8 sinh2g(t), 

which can be reduced to a particular case of the Painleve III equation. The connection 
with Fredholm determinants was discovered by McCoy, Tracy and Wu ||, and in 
the same paper the asymptotics as t — > of these solutions q(t) were determined. 
(Note that all asymptotics as t — > oo are trivial.) The asymptotics as t — > of 
det (J - X 2 K$) = det (I - A^ ) det (I + XK ) were determined in [[|. (See also 
0, where the asymptotics were found for a family of kernels including this one as a 
special case.) The asymptotics of det (/ — XK ) itself were stated without proof in 



A class of periodic cylindrical Toda equations arises in thermodynamic Bethe 
Ansatz considerations |§. There the additional constraint q~k-i = —qk is imposed. 
The solutions ( |1.2|) satisfy this constraint as long as the coefficients satisfy c w = 
—uj 3 c w -i. (This follows from the fact that det (I — K k ) = det (/ — K- k -2) i n this case, 
which is proved by applying the change of variable u — > u^ 1 .) The case n = 3 of this 
gives the cylindrical Bullough-Dodd equation (q = q% now) 

g /, (t)+t-V(*)=4(e 2 «W-e-«W), 

which can be reduced to another special case of Painleve III. Asymptotics of a class 
of solutions to Pjjj including this one were announced in || . 

This paper is devoted to the determination of the asymptotics of the quanti- 
ties det (J — XKk) in the general case, under the condition stated below. (In the 
final sections we shall compare our results in the cases n = 2 and 3 with those 
cited above.) We write K for Kq and consider at first only the asymptotics of 
det (/ — XK). This is no loss of generality since K k is obtained from K upon re- 
placing the coefficients by u k c w . The problem reduces to the asymptotics of 
Jq 00 R(u, u; X) du as t — > 0, where R(u,v;X) is the resolvent kernel of K, the ker- 
nel of K (I — XK)^ 1 . Using operator techniques, we show that R(u,u;X) is well- 
approximated on [1, oo] by the corresponding function when the exponentials in ( |1.3j ) 
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are replaced by e - ^ 1- ^" and on [0, 1] by the corresponding function when the expo- 
nentials are replaced by e~^ 1_w > u . (Actually the kernels have to be modified first 
by multiplying by factors (u/vf with ft depending on A.) We shall show that after 
these replacements we obtain operators which can be transformed into Wiener-Hopf 
operators, whose resolvent kernels have explicit integral representations. By these 
means the problem becomes that of determining the asymptotics of certain integrals. 
This is achieved by contour-shifting, and we find in the end that as t — > 

det(I-XK) ~b(~Y (1.4) 



with a and b constants given explicitly in terms of certain zeros of the function 

h(s) := sin ns — A n c u (— cj) s_1 . 

These are the values at A of those zeros which equal 1, • - • , n when A = 0. 

To state the result precisely, we denote by = ctfc(A) {k e Z) the zeros of this 
function indexed so that ctfc(0) = k. The zeros depend analytically on A as long as 
they are unequal, and when A = they are the integers. We derive the asymptotics 
( |1.4|) under the assumption that there is a path in the complex plane C running from 
to A such that everywhere on the path 

3?a <3?ai, $ta < 1, 3?«i>0, (1.5) 

and no zero lies in the strip 3?«o < < 3?«i. With this assumption the constants 
a and b are given by the formulas 

1^ 2 (n + l)(2n + l) 

a = - E a - ^ 

where a and a' run over the set {ai(A),- --,a n (A)} and G denotes the Barnes G- 
function 0. 

From these formulas we obtain the asymptotics of the solutions (|1.2|) . The re- 
quirement now is that everywhere on a path from to A we have for all k 

Kctfc < 9ft Qfc+i, k — 1 < 3? < k + 1. 

If this holds then 

q k {t) = Alog(-)+logS + o(l), 
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where for k — 1, • • • , n the constants A and .B are given by 



A = 2{a k -k), B= J] 



i<j<k 




k<j<n 



n 



r( 



71 



+ 1) 



n 



n 



and for other values of k are given by periodicity. 

As for the correct range of validity of the formulas, we conjecture that it is enough 
that 3ft «fc < Sftafc+i for all k for some path from to A, and that the extra condition 
k — l<dtotk<k + lis automatically satisfied then. In the cases which we consider 
in detail this is so and we obtain the correct range of validity. Another way of stating 
the condition is as follows. Define A to be the complement of 



Then the region of validity should be the connected component of A containing A = 0. 
The region for which we prove the result is the largest connected subset of this set in 
which the extra condition holds. 

Remark. It is shown in || that the more general class of kernels 



gives a solution to the the cylindrical Toda equations by the same formulas. Here p 
can be any finite complex measure supported on a compact subset fl of 



This assures that the operator is trace class. In case Q is the set of nth roots of unity 
other than 1 the condition is satisfied and the solution will clearly be n-periodic. 
We shall actually do everything in this more general case and we find asymptotic 
formulas of the form (|1.4j ) for the corresponding determinants, with the constants a 
and b being given by integral formulas involving the function h(s) now defined by 



In the periodic case h(s) is itself periodic and the integrals are expressible in terms of 
its zeros in a strip of width n. This is why the result here is so explicit. The require- 
ments on the afc stated above now refer to this function, and we assume throughout 
that they are satisfied. 



{A : 3?«fc = 3?«fc + i for some k). 



(1.6) 




{uEC: ftu; < 1, Mu- 1 < 1}. 




(1.7) 
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2. The approximating operators 



Recall that K is the operator on L 2 (R + ) with kernel 

r e -t[(i-^)«+(i-w _1 )n" 1 ] 
K(u, v) — j dp(u) 



—LOU + V 

We denote by R(u, v; A) the resolvent kernel of K, the kernel of R\ := K (I — A K) . 
It is well-known that 

d f°° 
— — log det (I — A K ) = / R(u,u;X)du, 
dX Jo 

the trace of the operator R\. Hence 

p\ poo 

log det (I- XK) = - / / R(u,u;/j,)dudp. (2.1) 

JO J 

In this section we are going to find a good approximation to the integral 
Jo°° R(u, u; A) du when A satisfies the condition stated in the Introduction. (After- 
wards we shall replace A by \i and integrate with respect to p over the path from to 
A throughout which ( |1.5| ) holds.) We begin with an observation. If we multiply the 
kernel K(u, v) above by (u/vf for any (3 then the resulting kernel still represents a 
bounded operator on Zv2(R + ) because of the decay of the exponential factor at and 
oo and, although the resolvent kernel changes, its value on the diagonal u = v does 
not. We are going to find approximations to the resolvent kernels for these modified 
operators, and precisely which (3 we take depends on A. 

Here is how we choose it. It follows from our main assumption that for each A 
there there exists s\ G (0, 1) such that the function h(s) given by (|1.7f) has no zeros 
on the line $fts = s\. (In fact the assumption guarantees that s\ can be chosen to 
vary continuously with A.) With this s\ we set (3 = | — s\. Notice that |/3| < |, a 
fact we shall need in order to apply our approximation argument. 

We write the kernel as 

KM = -f / dp(u) (2.2) 

K v J Jn —uju + v 

and denote the operator itself by K t . We do not display the dependence on /3, which 
is fixed for now, but use the subscript t to help the reader distinguish those operators 
that depend on t from those that don't. Both kinds will arise; the former have the 
subscript t and the latter will not. 

The exponential in (U) is the product e -K^)u e -t{i~^) u -\ For m > 1 the 
second factor is uniformly close to 1 when t is small while for u < 1 the first factor is 
uniformly close to 1. This suggests that the operators Kf with kernels 

Kt(u,v) := (-Y / —dp(u), Kf(u,v) := - / & - —dp(u,), 

x v / Jn — uju + v y v J Jn —uju + v 
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should in some sense approximate Kf on u > 1, u < 1, respectively, and therefore 
the resolvent kernels of these operators should approximate the resolvent kernel of K t 
on these intervals. We shall show that this is so, and that if Rf(u,v; A) denote the 
resolvent kernels of Kf, also on R + , then 

/ R(u,u; A) du = / Rt(u, u; A) du+ o(l), 
jo ' Jo 

(2.3) 

/OO /"OO 
R(u,u;X)du — J Rf(u,u;X)du + o(l) 

as t — > 0. We denote by P + multiplication by the characteristic function of (1, oo) and 
by P~ multiplication by the characteristic function of (0, 1). We shall use the notation 
Oi(A) to denote any family of operators whose trace norms are at most |A| times a 
function of t which is o(l) as t —>■ 0. (The subscript 1 refers to the trace norm. We 
shall also use the obvious notation oi(l) later on.) The main approximation statement 
will be 

P ± (I -XK t y 1 P ± = P ± (I-XKf)- 1 P ± + o 1 (X). (2.4) 
Relations (|2.3|) with A = 1 follow from this since it may be rewritten 

P ± {(I - A Kt)' 1 -I]P ± = P ± [(/ -XKf)- 1 -I]P ± + 0l (A), 

and if we take the trace of both sides and divide by A we obtain 



Here is an outline of the proof of (2A)- We use the matrix representations of our 
operators corresponding to the decomposition of L 2 (R + ) as the direct sum of the 
spaces L 2 (0, 1) and L 2 (l, oo). Thus (the equal sign meaning "has matrix representa- 
tion") 



I-XP-KtP- -XP'KfP 



I-XK, 



-A P + K t p- I -X P + K t P^ 

Because the nondiagonal corners of the matrix have the mutually orthogonal pro- 
jections P ± occurring as they do we will be able, with error c»i(A), to replace the 
operator K t appearing there by the operator Kq obtained from it by setting t = 0. 
Thus Kq has kernel 

K (u, v) = (-Y f ^— dp(u). (2.5) 

K v ' Jn —uju + v 

(Note the lack of consistency with the notation K^ in the introduction; this should 
cause no confusion.) If the diagonal entries I — X P ± K t P ± are invertible, we can write 
the resulting matrix as the product of 

I-XP-K t p- \ 

I-XP + K t P+ J 
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on the left and 

/ / -A (I -XP-K t p-)- l p-K P^ 

K -X(I -\P + KtP + )~ 1 P + K Q p- I 

on the right. Next, because of our assumption we shall be able to show that the opera- 
tors I—X P ± K t P ± are uniformly invertible for small t (i.e., the operator norms of their 
inverses are bounded) and that their inverses converge strongly to (I — XP KqP ) _1 
as t — > 0. (Recall that A t is said to converge strongly to A as t — > if A t f — > Af 
for all / in the underlying space.) This is actually the crux of the proof. After that 
it follows, because P ± KqP t are trace class, that with error 01(A) we can replace 
(/ — AP ± i^ t P ± ) _1 by (/ — A P^KqP^) -1 in the nondiagonal entries. Thus, if we 
define 



M :-- 



I 



V -A (/ - A P + KtP + Y 1 P + K P' 
we will have shown 

/ / - A P~K t P- 



-A (/ - A p-KoP-^p-KoP" 
I 



I-XK t 



V 



o 




M + o x {X). 



From this, using the uniform invertibility of the J — A P ± K t P ± again and the invert- 
ibility of the constant matrix M. (which we have to prove) we deduce 



(I - XK t 



-i 



M 



-i 








(J - A P + K t P~ 



+ o 1 {X). 



Now (here is the trick), applying an analogous procedure to the operator family 
/ — A Kf gives 



(I-XK+)- 1 = M- 1 








(/- XP+KtP+y 1 



+ 0i (A). 



It is clear that the lower-right entries of the two matrix products are the same, and 
this is exactly the statement 

P + (I-X Kt)- 1 P + = P+(I-X K+y 1 P + + oi (A), 

which is half of ( |2.4| ). The other half is obtained similarly. 
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Carrying out the details of the proof of this will require, first, some general facts 
about families of operators on a Hilbert space. 



Fact 1. If At converges strongly to an invertible operator A as t — > and if the A t 
are uniformly invertible then A^ 1 also converges strongly to A^ 1 . 

This follows from the assumptions and the identity A^ 1 — A" 1 = A^iyA — At) A -1 . 
The next fact says that strong convergence can sometimes be converted into trace 
norm convergence. 

Fact 2. If A t —>■ A strongly and B t — > B in trace norm then A t B t — ► AB in trace 
norm. 

This is a variant of Proposition 2.1 of ||. There the families of operators depended 
on a parameter n G Z + rather than t G R + , a matter of no importance since we may 
consider general sequences t n — > 0. Also, instead of a sequence of operators converging 
in trace norm to B there was the single trace class operator B. The apparently more 
general result follows trivially from this special case. 

Fact 3. Suppose A t and A are as in Fact 1 and that B t are trace class operators 
converging in trace norm to B. Assume also that A + B is invertible. Then the A t + B t 
are uniformly invertible for sufficiently small t, and if B t = oi(l) then (A t + B t )~ l = 
At 1 + o 1 (l). 

Proof. Write A t + B t = A t (I + A^ 1 B t ). By Fact 1 A^ 1 -> A' 1 strongly, and so 
by Fact 2 (with A t repaced by A^ ) we deduce I + A^ x B t = I + A^B + oi(l) = 
(/ + oi(l))(7 + since clearly (J + 0!(1)) _1 = J + Oi(l). Both statements now 

follow. 

In our derivation of ( |2.4| ) we have to know that the / — P ± KfP ± are uniformly 
invertible for small t. We shall deduce this from known facts about uniform invert- 
ibility of truncated Wiener-Hopf operators, which we now describe. The proofs can 
be found in [|J. 

The Wiener-Hopf operator W associated with a function k G L\ (R) is the operator 
on L2(R + ) with kernel k(x — y). Introduce the Fourier transform of k, 




This is a continuous function on R tending to as £ — > ±oo. A necessary and 
sufficient condition that I — W be invertible is that 1 — ^ for all £, and 



oo 



arg(l-A;(0) 



0. 



— oo 
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The truncated Wiener-Hopf operators are the operators P a WP a where P a denotes 
multiplication by the characteristic function of (0,a). Clearly these operators con- 
verge strongly to W as a — ■> oo. The important fact is that if W is invertible, in other 
words if the conditions on k stated above hold, then the operators / — P a WP a are 
uniformly invertible for sufficiently large a. We mention also that the operator with 
kernel k(x — y) on the whole line R is invertible if and only if the first condtion alone 
is satisfied, that 1 — £;(£) 7^ for all £. 

From this we can deduce information about kernels k(u, v) which are homogeneous 
of degree —1 since the variable change u = e x transforms this kernel into a convolution 
kernel. More precisely, denote by U the unitary operator from L 2 (l, 00) to L 2 (0, 00) 
given by Uf(x) = e x ^ 2 f(e x ). Then if T denotes the operator on L 2 (l, 00) with kernel 
k(u,v), the operator UTU' 1 is the operator on L 2 (0, 00) with kernel 

e x/2 e v/2 k(e x , e y ) = e (y ~ x)/2 k(l, e y ~ x ), 

where we used the homogeneity of k(u,v). Notice also that if Pf denotes multi- 
plication by the characteristic function of then UPfU~ l is the projection 
operator Pi og t-i of the last paragraph. After making an obvious change of variable in 
computing the Fourier transform of e~ x l 2 k(l, e~ x ) we deduce 

Fact 4. Assume that J °° v ~ 1//2 \k(l, v)\dv < 00, denote by T the operator with kernel 
k(u, v) on L 2 (l, 00) with kernel k(u,v), and by M(s) the Mellin transform 

M(s) := / V s ' 1 k(l,v)dv. 
Jo 

Then a necessary and sufficient condition that / — T be invertible is that 

l-M(s) ^Ofor^s = i arg (1 - M(s)) = 0. (2.6) 

- 2 -ioc 

If this holds then the operators / — P^TP^ are uniformly invertible for sufficiently 
small t. 

Remark. If we use the variable change u = e~ x instead of u = e x then our operator 
T acts on L 2 (0, 1) and we find (again using homogeneity) that the condition for in- 
vertibility of I — T is exactly the same as before, and that if Pf denotes multiplication 
by the characteristic function of (t, 1) then this condition implies the uniform invert- 
ibility of I — PfTPf for sufficiently small t. Also (tranferring to this context the last 
sentence of the discussion of Wiener-Hopf operators), the same condition implies the 
invertibility of the operator / — T on L 2 (0, 00). 

We apply this to the kernel \K , which is homogeneous of degree —1. The relevant 
Mellin transform is found to be 



M(s) = A — - — / {-uy-^dpiu). 
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so that 1 — M (s) = h(s — j3), where h is given by (|1.7|). If we recall that (3 = \ — s\ we 
see that the conditions ( |2.6|) are met, the first immediately from the definition of s\ 
and the second because the index (the variation of the argument , which is necessarily 
an integer) is a continuous function of A, locally constant and clearly equal to when 
A = 0. So we know that the operators / — A Kq and / — A P ± KqP ± are invertible and 
the operator families I — X P^K P^ are uniformly invertible for small t. 

We introduce one last piece of notation. We denote by the operators on 
L 2 (R + ) with kernels 

K + (u,v) = P.) / ^—dp(co), K-(u,v) = P.)* / G - —dp( U ). 

K v J Jn—uiu + v K v J Jn —uju + v 

Notice that rescaling K ± (u ) v) under the variable change u — > t ±l u gives K^(u,v). 

Next we derive various trace class properties of our operators which will be needed. 
For these shall use an estimate for the trace norm of an operator on L2(R + ) with 
kernel of the form 

qi (u,u) q2 (u,,v) dp{u;) (27) 



— urn + v 

which is a special case of the sublemma in the appendix of ||. We denote by 0(s) any 
positive function on R + , by $ its Laplace transform, and by \1/ the Laplace transform 
of Then there is a constant m depending on only on Q such that the trace 

norm of the operator with kernel given above is at most m~ l times the square root of 

/•oo r- 

qi(u,u)\ 2 $(mu) dp{uj) du ■ / / ^(w, u)\ 2 ^(mu) dp(ui) du. (2.8) 

J o Jn 

From this will follow our first lemma. In the proof we denote by X + the charac- 
teristic function of (1, oo) and by X~ the characteristic function of (0,1), so P^ is 
multiplication by X ± . 

Lemma 1. The operators (independent of t) 

P + K P-, P+K+, P~(K+-K ) 

are trace class. The operators (depending on t) 

P + (K t -K )P~, P-(K+-K ), P + (K t -Kt) 

are Oi(l) as t — > 0. The statements also hold if all superscripts + and — are inter- 
changed. 
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Proof. All the operators have kernel of the form ( |2.7D . We list the operators below, 
together with the corresponding functions qi(u,u) and q 2 (uj,u). 



Operator qi(u,u) q 2 (uj,u) 

P + K P- X {1>oo) {u)u p X {Q ,i){u)u- p 

P+K+ X(i,oo)(V)e~ (1 ~ a;)u u /3 u~P 

P-(K+ - K ) Xmi)(u) (e-^^ - 1) vP u 



P-{Kt - K ) X m (u) (e-^-^ u - 1) vP u'f 3 
P + (K t -Kt) Xji^He-'^'^e-^" 1 !''" 1 - l)u p 

For each of these operators we take two numbers p,q £ (—1, 1) and define 4>(s) = s 
for s < 1 and 0(s) = s 9 for s > 1. We easily see that 

OIu-p- 1 ) for u > 1 f O^- 1 ) for u > 1 

4>f.s) { V(s) = \ 

(Din-i- 1 ) ior u < 1, Ofw 9 " 1 ) for u < 1. 



For each operator one can find p and g such that both integrals in ( |2.8| ) are finite, and 
any integral depending on t is o(l) as t —>■ 0. In fact, as the reader can check, we may 
take for all the operators any q £ (2/3, 2/3 + 2), for the first and fourth operators any 
p > 2(3 and for the other four any p £ (2/3 — 2, 2/3). This takes care of the six displayed 
operators. For the other six we use the fact that the substitutions u — > u~ x , c<j — > uo~ l 
yield operators of the same form with (3 replaced with — j3 and with the superscripts 
interchanged. 

As a preliminary to the next lemma we show that certain modified Laplace tran- 
forms are bounded operators on L 2 (R + ). 

Lemma 2. The integral operator on L2(R + ) with kernel (uv)~@ e~ uv is bounded if 

(3<i 

Proof. The mapping f(v) — > v^fiv^ 1 ) is unitary, so we may replace the kernel 
by u~^v^~ x e~ u l v . Under the unitary mapping f(u) — > e x ^ 2 f(e x ) this becomes the 
kernel e^~ l3 ^ x ~ v 'e~ eX v on (— oo, oo). Thus the operator becomes convolution by the 
Li function k(x) = e^~^ x e~ eX and so is bounded. 

Lemma 3. The operators I — A K are invertible. 

Proof. We consider K + , which we can write as P~K + P~(K + — K ) + P + K + . By 
the first part of Lemma 1 the second summands are trace class, therefore certainly 
compact. Our assumption implies that / — \P~K is invertible. Hence / — \K + 
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is the sum of an invertible operator and a compact operator, and so it follows from 
general theory that it will be invertible if is not an eigenvalue. In other words it 
suffices to prove that \ K + f = f for / G L 2 (R + ) implies / = 0. 
For any uj G C\R + and any / G L 2 (R + ) we have for x > 



e- xu u-?du I (-Y 6 (1 — f(v)dv 
o Jo K V J —uju + v 



dy 



lo x + 1 -u(y+ 1) Jo w 
This can be seen for 3ft uj < by using the integral reresentation 

1 -yy— 



—UJU + V 



in the integral on the left above and interchanging the order of integration. The 
identity follows for all uo G Q since both sides are analytic functions of uj in this 
domain. Now suppose that XK + f = f. Then if we integrate both sides of the 
identity with respect to dp(u) and multiply by A the left side becomes the Laplace 
transform of m _/3 /(m), which we denote by g(x), and the the right side becomes 

Jo Jnx + l-uo(y + l) yKyj y 

Using the fact y~^g{y) G L 2 , which we know by the previous lemma, we see that the 
integral is a bounded function of x. If we recall the definition ( |2.5| ) then we see that 
the identity becomes 



poo ,x A- 1 \ 3 

X (——) K (y + l,x + l)g(y)dy (x > 0), 
Jo K y + l J 



or 

x- p g(x - 1) = A j~~ K (y, x) g(y - 1) y^ dy (x > 1). 

Now we know that x _/3 g(x) is in L 2 (0, oo) and that g{x) is bounded. It follows that 
x~^ g{x — 1) belongs to L 2 (l,oo). The right side above is the operator with kernel 
K Q (y, x) acting on this function. Thus (if / ^ 0) the operator A P + K' Q P + has 1 as an 
eigenvalue, where ' denotes transpose, so I — A P + K' Q P + is not invertible. But this 
implies I — \P + K P + is not invertible, whereas we know that it is. This contradiction 
establishes the lemma. 

Lemma 4. The operators I — X P ± K^ 1 P ± are uniformly invertible for sufficiently 
small t. 
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Proof. We consider I — X P + K^P + and for this it is enough to show that the 

1 — \P + Kt~ are uniformly invertible. The kernel of P + K^P + is X(i t00 )(u) Kf(u, v ) 
and the substitution u — > t~ l u allows us to consider instead the operator with kernel 
X(t,oo)(u) K + (u,v). We write this (not displaying the variables u and v) as 

X {t ,i) K + X ( i i0o) K + + X (0 ,i) (K+ - K ) + X m {K - K + ). 

Recalling the definitions of our various projection operators we see that the first kernel 
corresponds to the operator PfK , and we know that the I — X PfK are uniformly 
invertible for sufficiently small t. The second and third summands correspond to the 
operators P + K + and P + (K + — Kq), which we know by Lemma 1 to be trace class. 
The last summand corresponds to the operator P~(Kq — K + ), which we know by 
Lemma 1 to be trace class, left-multiplied by multiplication by X(o,t) , which converges 
strongly to 0. An application of Fact 2 shows that this last operator is Oi(l). The 
strong limit of the sum of the four operators is, of course, K + and we know by Lamma 

2 that / — \ K + is invertible. Hence we can apply Fact 3 to deduce the the result. 



We can now fill in the details of the proof of ( |2.4| ) outlined earlier. Thus we begin 
with the representation 



I-XP-K t P- -XP-K t P 



I-XK f . 



-A P + K t p- I -X P + K t P^ 

Applying Lemma 1 to the nondiagonal entries we deduce 

fl-XP'KtP- -XP-K P + \ 
I-XK t =\ +oi(A). 
V -XP + K P- I-XP + KtP + ) 

Lemma 3 tells us in particular that the diagonal entries of this matrix are invertible 
for small t so we may factor out 

I-XP-K t p- \ 

I-XP + K t P + ) 

on the left, leaving 

/ / -A (I - XP-K^p-^p-KoP 4 

K -X (I - X P+KfP+^P+KoP- I 

on the right. Next we combine the uniform invertibihty of the / — A P ± KfP ± proved 
in Lemma 3 with Fact 1 to deduce that the inverses of these operators converge 
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strongly to (I - XP ± K^P ± )~ 1 . S ince P KqP^ are trace class, by Lemma 1, we 
deduce by Fact 2 that the matrix above is M. +oi(A) where M. is the matrix obtained 
by replacing Kf by K . Thus 



/ I - X p-KtP- 



I-XK t 







V 



o 



I-XP + K t P + J 



M + 0l {X). 



Now we have to know that Ai is invertible, and we see this as follows. If, instead 
of the operator I — XK t which depends on t, we had started with the operator / — Kq 
then we would have obtained the exact representation 



/ / — XP~KnP~ 



I-XKn 







\ 



V 



o 



M. 



I-XP + K P + J 



Since both I — X Kq and the matrix on the left are invertible, by our assumption, we 
deduce that M. is invertible. 

Next we go through a similar process starting with the operator I — X Kf rather 
than 1—XK. Using the fact that P~Kf = P~ K +oi(l), which we know by Lemma 1, 
we obtain in this case 



/ I-XP-K n p- 



I-XKt 



\ 











/ - A P + K t P + 



M + 0l (X). 



From these matrix representations and the facts that M. and / — XP~K P~ are 
invertible and / — P ± KfP ± uniformly invertible we deduce, using Fact 3 with Oi(l) 
replaced by Oi(A), 



(I - XK t 



M 



-i 



(/ - XP-KtP-y 1 

(I -XP + K t P + y l 

(i - xp-KoP-y 1 

o (i - xp+KtP+y 1 J 

Comparing lower-right entries of the matrices gives 

P + (I - K)- 1 P + = P + (I- KtY 1 P + + oi(A), 



(I-XK t 



M 



-i 



+ oi (A), 



+ oi(A). 



which is half of (2~4). The other half is obtained similarly. 
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Remark. To apply ( |2.3| ) to ( [2.1D we need something extra, e.g., that ( p.3| ) holds 
uniformly for these A. With a little care our argument gives this also, but we spare 
the reader the details. 

3. The resolvents of 

We are going to find integral representations for the integrals on the right side of 
(P- 3D, and we consider Rf~(u, u; A) du first. The substitution u — > u/t shows that 
this equals f£° R + (u, u; A) du where R + (u, v; A) is the resolvent kernel of the operator 
K + . For this we require only that Q be a compact subset of 

{lo G C : ftu; < 1, io G'R+j, (3.1) 

since the term 1 — u^ 1 does not appear in the exponent in the kernel of K + . The 
derivation will involve an initial step which is valid only when Q is contained in the 
left half-plane so we assume this to begin with. We shall also assume that A is so 
small that h(s) ^ for 3?s = ~, so that with the notation of the last section we may 
take S\ — |, (3 = 0. Eventually these two assumptions will be removed by an analytic 
continuation argument. 

Because (3 = the kernel of K + is 

K + (u,v) = / —dp(u). (3.2) 

If we set 

A(u,x):= e-^-^e^dpito), B(x, u) := e~ ux , 
Jn 

then _ 

POO 

K + (u 1 v)= / A(u, x) B{x, v) dx. 



Lemma 2 of the preceding section tells us that B(u,x) is the kernel of a bounded 
operator from L2(R + ) to L2(R + ) and, with our assumption on Q, that the same is 
true of A(x,u). The above shows that K + = AB, and the operator BA has kernel 



B(x, u) A(u, y) du 



dp(uj) f dp(uj) 



n x — uy + 1 — to Jn (x + 1) — uj{y + 1) 



We use the general fact AB(I - XAB)- 1 = A(I - \BA)- 1 B to deduce that R + (u } u) 
is given by an inner product, 



R+{u, u) = ((/ - \BA)- l B{ -,u), A{u, ■ )) . (3.3) 
We begin by computing 

/:=(/-A J BA)- 1 J B(-,n). 
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Thus we want to solve 

m -a/7 - ( f(v) d y = e ~ ux ( x * o), 

Jo ./n (x + 1) — u(y + 1) 

or 

f(x - 1) - A r / -^H /(y - 1) dy = e"^- 1 ) (x > 1). 

The substition x — > e x brings this to the form of a Wiener-Hopf equation, so we can 
use the factorization method to find the solution. 

We begin by decreeing that the last indentity holds for all x > 0, in other words 
we define / on (—1,0) by the identity. Then we define 

/•oo rl 

F_(s):= X s ' 1 fix - 1) dx, F+{s) := / X s " 1 f(x - 1) dx. 
Ji Jo 

These belong to the spaces Hardy spaces H 2 (JRs < |), H 2 (?Rs > |), respectively. We 
take Mellin transforms of both sides of the equation, and find that for 9ft s = | 

F_(s) + F+(s) - ( (-uj) s - l dp(uj)F_(s) = e u u- s T(s). 

sin its Jn 

(The exponential in the integral is made definite by taking | arg(— u)\ < n.) We write 
this as 

H(s) F^s) + F + (s) = e u u- s T(s), (3.4) 

where 

H(s) : = = 1 - I (-uy- 1 dp{cu). (3.5) 



sm its sin its Jn 

This function is bounded and analytic in each vertical strip of the complex s-plane, 
away from the zeros of sin its, H(s) — 1^0 exponentially as — > ±oo and 



arg H(s) 
Thus there is be a representation 



His) H - (S) 



H + (s) 



where H_(s) ±1 are bounded and analytic in 9ft s < \ + 5 for some 5 > and H + (s) ±1 
are bounded and analytic in 9ft s > | — 5. We multiply (|3.4|) by H + (s) and use the 
decomposition F = F_ + F + of an arbitrary function in L 2 (3fts = |) into boundary 
functions of functions in H 2 ($ls < |) and H 2 (^fts > |) to write the result as 

if_( S ) F_(a) - e" ( M - s r(a) = -H + {s) F + (s) + e M (V s T(s) H+{s)) 
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The two sides are boundary functions of functions in if 2 (3i!s < \) and H 2 (Jfts > |), 
respectively, so they both vanish. This gives the representation 



F4s) = jr ^{u- s T(s)H + (s))_. (3.6) 

Now (see (|3~^)) we have to multiply f(x) by A(u,x) and integrate with respect to 
x over (0, oo). This is 

dp(u) / f(x) e - (1 - w)M e" ux dx = e~ u / dp{u) / f(x - 1) X {1>oo) {x) e wux dx. 



The Mellin transform of f(x — 1) X(i,ca){ x ) equals F-(s) and the Mellin transform of 
6 ujux e q Ua i s uj U ^- s r(s), so Parseval's formula for Mellin transforms shows that the 
above equals 

e~ u [ I (-u;uy- 1 dp(u;)V(l-s)F4s)^ 
J Jn 2m 

the outer integration taken over 9ft s = |. (All vertical integrals are taken in the 
direction from — ioo to zoo.) Next we recall (|3.6j ) and use the integral representation 
of the operator G — > G_ to write the above as 

, .r(l-s) rfs f u~ s ' T{s')H + {s') ds' 



ln y rv ' H_(s) 2m' J s' — s 2m' 

the inner integral taken over 3ft s 1 = ~ + 5. Alternatively, this may be written 
_ T(l-s) rfs /•M- s '- 1 r(s' + s)F + (s / + s) ds' 



n 



Hjs) 2m J s' 2m' 



where now the inner integral is taken over 3ft s' = 6. The integrands of these integrals 
vanish exponentially at infinity, and u occurs to the power —s' — 1, which has real 
part —5 — 1. Thus we may integrate with respect to u from t to oo under the integral 
signs and deduce that 



R + {u,u;X)du= / {-uj) s - 1 dp{u) 



. s _ x , , . T(l-s) ds ft- 3 ' 'r(s' + s)H+{s' + s) ds' 



n ' rv ' H_(s) 2m J s' 2 2m 

It follows from ( |3.5| ), and the gamma function representation of the last factor there, 
that 

' (-uy- 1 dp{u) r(i - s ) = a- 1 (i - h(s)) r{ s )-\ 



n 



Thus we have shown (reverting to the resolvent Rf) that 

POO 

A J R^(u,u;X)du 
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v #_(s) -if + (s) )r(s) — 



2iri Js-ioo s' 2 2%i 

This was proved if A is sufficiently small and if Q lies in the left half-plane. Let us 
remove the latter condition first. For any rj define the measure p v by p v {E) = p(E—rj). 
This has support Q + rj. For all rj in a neighborhood in C of [—1, 0] the set Q + rj is 
contained in the region (3T). For rj near —1 the set will also lie in the left half-plane. 
If A is small enough the condition on the zeros will be satisfied for the measures p v 
for all r\ in a neighborhood of [—1, 0]. For such A we know that the above formula 
holds for rj in a neighborhood of —1. But both sides are analytic functions of rj in our 
neighborhood. Thus the formula must hold for rj = also, which is what we wanted 
to show. 

To remove the condition that A be small we must modify the formula to read 



/oo 
Rf(u, u\ A) du 

l ds r s +^ t- s 'T(s' + s)H + (s' + s) ds' 



(II (.-) ' II. (s) ; )!•(*, , 

s\- 700 Z7TZ J 5—ioo S Z7T1 

(3.7) 

where S\ G (0, 1) is as in the previous section, a continuously varying function of A 
such that 1 — H(s) is nonzero on the line 3?s = s\. Both sides of ( |3.7|) are analytic 
functions of A for A in a neighborhood of our path, they agree near A = 0, so they 
agree everywhere on the path. 

As for R~(u, u\ A) on (0, 1) the change of variable u — > u^ 1 transforms the kernel 
K~(u, v) into 

„ p —t(l—U!~ 1 )u „ — t(l— Ul)u 

K + (u,v) := / —dp{u) = / — (-a;)^- 1 ). 

Jn —ujv + u Jn —iou + v 

Therefore 

R-{u, u; A) = u~ 2 R + (u-\ u" 1 ; A), 
where R + is the resolvent kernel for K + , and so 



pi roc „ 

/ R~(u, u; A) du = / R + (u,u; A) du. 
Jo Ji 



It is easy to see that replacing K + by K + replaces H(s) by H(l — s) and so the 
integral is equal to (|3.7|) with H(s) replaced by H(l — s) and s\ replaced by 1 — s\. 



Now that we have these explicit representations it is obvious what we do: in the 
inner integral in (37?) we move the line of s'-integration from 3? s' = 5 to 3? s' = —5. 
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We can do this if 5 is small enough. The residue at the double pole at s' — 
contributes 

- I (H(s)- l -l)^\ogt+ [ (H^sr'-H+isy^-^ir&H+is))'^ (3.8) 
J lixi J 1 [s) Lix% 

the integrations taken over 3?s = s\, and the error term is 0{t s ). We do the same 
with H(s) replaced by H{1 — s) and add, and so we have obtained the asymptotics 
of A / °° R(u, u\ A) du. 

4. Asymptotics of det (/ — K) — The periodic case 

The formula for H(s) is now 



sin Tis V 



uj running over the nth roots of unity other than 1. If in our sums we set — u = 
e v( 2 -?~ n ) ; (j = 1, • • • ,n — 1) then | arg (— u>)\ = \f t (2j — n)\ < tt as required. If we also 
set z = e^ s then the above may be written as 



n-l 

2j-n 



sin us H(s) = sin its + Xrr V] cu 1 c„, z 



or 

n-l 



2i sinvrs H(s) = z~ n z 2n - 1 + 2i\n £ - 



i=i 

Recall that smirs H(s) = h(s). The expression in brackets above is a polynomial of 
degree n in z 1 and its zeros are the quantities by e~^ a , where a runs through the 
zeros ojfc = ctfc(A) (/c = 1, • • • , n). With this notation the right side above is equal to 

z -n JJ(^ - e^ a ) = H(ze~ v ° _ 2 -i e >) J] e S«. (4.1) 

a a a 

Here and below the index a runs over the set {aci, • • • , «„}. The last product, a square 
root of the product of the roots of the polynomial, equals ±1 or ±i. This product 
equals e^ 1H hn - ) = i n+1 for A = and so for all A. Recalling that z = e^ s we see 
that we have obtained the representation 

H( s ) = K -^>-± TT sin l( s _ Q,). 

sm7rs „ n 



We now evaluate the integrals in ( |3.8| ). For A sufficiently small again, the will 
all lie in the strip | < 3?s < n + ~. We may assume this since the usual analytic 
continuation will give the general case. 
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To evaluate the first integral in ( |3.8| ) we consider 



ds 
2m 



taken over the infinite rectangle which is the contour running from n + | — ioo to 
n + \ + ioo and then from \ + ioo to \ — ioo. On the one hand this equals n times 
the first integral in ( |3.8| ), and on the other hand it equals the sum of the residues at 
the poles between the two lines. Thus we have shown that 

(^)- 1 -l)^ = ^E«^'(«)- 1 - (4-2) 



For the second integral in ( |3.8| ) we have to write down the explicit expression for 
the factors H±(s). These are given by 



r(«) ~ w n„r(- 



s — ce i 
n 



It is readily verified that H(s) = H_(s)/ H + (s) and that if_(s) ±:L and if + (s) ±1 are 
are bounded and analytic in 3ft s < \ + 5 and 3ft s > | — 5, respectively, for small A. 
Thus, they are the correct factors. The second integral in (|3.8|) may be written 



/ 



_! ^(T(s)H + (s))' ds 



His)' 1 -! 



T(s)H + (s) 2m 
and by above expression for H + (s) this equals 

/I F'f s ~ a ' _|_ 1 

(^(s)- 1 - 1) [- y ; " , ; + i og n 



(is 
2^1' 



where in the sum a' also runs over the set {cei, • • • , a n }. The contribution of the term 
logn is exactly logn times (|4.2| ). We evaluate the rest of this integral we use the 
characteristic property of the Barnes G-function, G(z + 1) = T(z)G(z). Putting z 
equal to + 1 and taking logarithmic derivatives gives 



We integrate 



G'(^ + 2) G'j^f + 1) _ T'j^f + 1) 
G(^ + 2) G(^ + l) r(^ + i)' 

n'l s-a' I i \ 

(fr(a)- 1 - 1) £ " 



G(^ + l) 

over the same infinite rectangle as before. (This is justified by the fact that H(s)^ 1 — 1 
vanishes exponentially at oo in vertical strips while G'(z)/G(z) grows like z log z.) 
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By the above relation the result is exactly the integral we want, and so computing 
residues gives the formula 

1 f , r , n 1 x r '(— + 1) ds 1 ^ G'(s=sL + i) 



[ (H(s)- 1 -1)V 1 l~ + ^ = £^ 

n7 1 lJ ; rf^si + i) 27ri n^Gf^ + l) 

a' \ n I ol,o! \ n > 

Thus, we have shown that 



A 

where 



/ R+{u, u, A) = a + (A) log (-) + 6+ (A) + 0(t s ), 



a+{\) = - l -Y.*H\a)-\ ^{X)^W %Zl 1 } H'{a)-\ (4.3) 

We must add to this A Jq 1 Rt{u, u, A) du which, as was mentioned earlier, is obtained 
by replacing H(s) by if (1 — s). The zeros of this function which lie near 1, • • • , n for 
small A are n — a + 1. Hence (|4.3|) is replaced by 

a" (A) = -Y(n-a+l) H'(a)-\ 
n ^ 

b- (A) = -- V G ^ + \ } iy'H" 1 = -- V G ^ + l ] H'(a')- 1 - 

Here we have used the periodicity of H and the fact £H(l — s) = —H'(l — s). Adding 
and using (|2.3|) , we see that 

X R(u, u; A) <iw = a(X) log ( - ) + 6(A) + 0(t ). (4.4) 
Jo K n J 

where a(A) = a+(A) + a" (A), 6(A) = 6+(A) + 6" (A). 

Now to obtain the asymptotics of log det (I — K) we must replace A by /x, multiply 
the above by —dfi/fi and integrate from to A. (Notice the factor A on the left side 
of Q4.4|) and recall the minus sign in (|2.1|) .) We obtain from (|3.5| ) that for a zero a (A) 
of H we have Xda/dX = H'la)' 1 . Thus the coefficient a(A) in ( f4.4| ) may be written 
(after replacing A by \x and thinking of a as a(fi)) 

Y (2a — n — 1) /i da/ da. 

n „ 



Multiplying by —d\ij\i and integrating gives 



1 



n 



Y(a 2 -(n+i)a) ;: =;e 



X 1 



m=o n 



a' 



fi=X 
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since, as we have already seen, a is independent of A. Similarly b(/i) may be written 

1 G^+l) , ... 

-E^EZTI)^ (a - a)/ ^ 

and multiplying by — and integrating gives 

E, „/Q! — a m=A 
logG + 1 
, n m=o 

a, a' 

If we recall that when \i — the zeros are 1, • ■ • , n we see that the formulas for the 
constants in (|1.4j) are the ones stated in the introduction. 

To obtain the asymptotics of the quif) we must consider det (I — XKk) instead of 
det (I — XK). This amounts to replacing the coefhients c u by u k c w , and this in turn 
amounts to replacing H(s) by H(s + k). The zeros of this function modulo n are 
ai(X) — k, • ■ ■ , a n (X) — k. But these are not the zeros which are to replace the a, a' 
in our formulas for a and b since they do not arise from from the zeros whose values 
are 1, • - • , n when A = 0. Rather, the replacements must be 

a k+1 (X) - k, ■ ■ ■ , a n (X) - k, ai+n- k(X), ■■■ ,a k + n- k(X), 

which are the zeros with this property. Thus, for the asymptotics of qk{t) we make 
these replacements in our formulas and the corresponding replacements with k — 1 
instead of k, subtract, and take logarithms. The result is found, after some computa- 
tion and the use of the functional equation for the G-function, to be the asymptotics 
stated in the introduction. 

5. Asymptotics of det (/ — K) — The nonperiodic case 

The coefficients a(A) = a + (A) + a~(A), 6(A) = b + (X) + b~(X) of the last section 
were in general given by integral formulas. They were 

a+(A) = - / (H(s)- 1 - 1) ^L, (5.1) 

J Z7T2 

6+(A)-a+(A)logn = J (H^ - H + (s)-') ^(T(s) H + (s))> (5.2) 

with the formulas for a~(A), and b~(X) obtained by replacing H(s) by if (1 — s). The 
integrations may be taken over 3?s = | if A is small enough and, as usual, this is 
no loss of generality. To find a and b we integrate a(/x) and &(/i), respectively, with 
respect to —dfi/fi over a path from to A. 
Write 

<p{*) :=-^ / (-uy-'dpiu), 
sin 7ts Jn 
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so that H(s) = 1 — X(p(s). Making this replacement in the integrand in ( |5.1| ) and 
integrating gives 

- f (~j ~~TT " 1) ^ = " f 4> = log (1 - «.)), 

so the contribution to the coefficient of log t is 

ds 



J log(l-A^( S )) 



27U 



over 3?s = ~. Replacing H(s) by — s) gives the same contribution since we may 
make the substitution s — > 1 — s. Therefore 

/ds 
log(l-A^(s)) — . (5.3) 

For general A the integration is to be on 3ft s = s\. 
Now we go to (|5.2|), which may be written 



, . , / r'(s) / 1 a ds r fH'(s) H'(s)\ ds .„ 

6 + (A) - a + (A) logra = / — )-f — — - 1 — : + / ( +v ( - +) '- ) — : . (5.4) 



r(s) v #(s) / 2m y V #-(s) #+(s) y 2m 

By a computation like the earlier one we see that the first integral becomes after the 
/i-integration 

Then we replace s by 1 — s, make the substitution s — > 1 — s, and add. We see that 
the contribution of the first integral in (|5.4j) equals 



/( 



r'(s) r'(i - s )n . „ , , „ ds 

+ r) log(l- \<p(s))— . (5.5) 



T(s) r(l - s) / ov rv " 2m 



Finally, we look at the second integral in ( |5.4j ), which equals 

f^M—=[ ^Mh(s)- 1 — = [ (log H + )>{ 8 ) H(sr — • 
J H_(s) 2m J H + (s) K J 2m J v 6 +; v ; v ; 2m 

Replacing H(s) by H(l — s) replaces H + (s) by 1/H-(1 — s), so after making the 
substitution s — ► 1 — s and adding we get 

/ (logif + iy_)'( S )iy( S )- 1 ^ : . (5.6) 

J 2717 

Recall that for 3ft s = | 

, tt , S 1 , rr / \ f log tf(s') ^' 
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where the integral is the Hilbert transform, a principal value integral over 3ft s' 
So 

i u ( \ u ( \ f lQ g H ( s> ) ds ' 
log H+(s)H-(s) = / — r. 

J s' — s m 

Since the Hilbert transform commutes with differentiation we get 

n ff ff vn 1 f H ' (s ' ] ds 
log H+HJ) (s) = — 

m J 

and so ( |5.6|) equals 



if(s') s'-s : 



2tt 2 y y H(s') 
The //-integration gives 

V( s ) Ts nWTi TTT^ 



1 — fKf(s') l-flif(s) fi JO (1 — pLlf(s')) (1 — (i(p(s)) 

if f\ 1 , 1 - V(s) 

= V ( s ) 7 a TT log 



y(s') — <£>(s) 1 — X(p(s') 
and so the contribution of the second integral in (|5.4j) is 



27T 2 y y v 9 ( s _ ^( s ) i _ ^v 9 ( s s' — s 



Thus 6 — a logn equals the sum of (|5.5|) and (|5.7|) . As usual, for general A the integrals 
are taken over 3ft s, s' = S\. 



Remark. The double integral (|5.7| ) is exactly the constant in the known asymp- 
totics for the determinants of the truncated Wiener-Hopf operators associated with ip 
(specifically, <^(| + is the Fourier transform of the convolving kernel), and (5.3) is 
(minus twice) the leading coefficient in the asymptotics. One can see by the argument 
of Section 2 how both these things arise and conclude also that Q5.7D equals det M.^ 1 . 
The extra ingredient here is therefore the integral (|5.5|) . 

6. The case n = 2 

In this case the only root is uj = — 1 and we may take c_i equal to 1 since it occurs 
only in the product Ac_i. Thus the kernel of Kq is 

-2t(«+«- 1 ) 



U + V 

and the equation (for either g*.) is 

q"(t)+r l q\t) = 8 sinh2g(t). (6.1) 



24 



We have in this case h(s) = sin7rs — tt\, the zeros are given by 

«o = — arcsin7rA = — \og(ni\ + \/T— 7r 2 A 2 ), ati = 1 — a , 

71 1X1 

and «fc+2 = afc + 2. The square root is that branch which is positive for A = 
and the logarithm that branch which is there. From this it is easy to see that 
the set fllDD consists of the rays (— oo, — l/7r] and [1/vr, oo) and A, the proposed 
region of validity of our formulas, is the complex plane cut along these rays. If we 
note that the function + y/1 — n 2 \ 2 map A onto the right half-plane, we see that 
|3?a!fc(A) — k\ < \ for all A G A and so the "extra" condition on the a is satisfied. The 
range of valididy is therefore all of A. Using the formulas stated in the introduction 
we find that det (J — A K ) ~ b (t/2) a with 



g%g% r(|) g(|; 



2 



G{\ + oo) G(f - a ) r(| - a ) G{\ + a ) G{\ - a ) ' 

For det (/ + XK ) we replace A by —A, which amounts to replacing a by — a . If 
we multiply the two results together we recover the asymptotics for det (/ — A 2 Kq) 
determined in and 0. 

For q we have the asymptotics Alog(t/2) + logB + o(l), where 

r(i — an) 

A = 2a , B 



r(| + a )' 

in agreement with [||. This is the solution of (|6.1|) which is asymptotic to — 2A-K" (4t) 
as t — * oo, where this i^o is the Bessel function. 

For A ^ A the asymptotics are different. For A > 1/tt e w ™ has an infinite sequence 
of zeros as t — > and for A < —1/tt it has an infinite sequence of poles; this follows 
from the fact that as t — > the spectrum of K fills up the interval [0, tt]. A heuristic 
derivation of the asymptotics for A on the cut is given in 0. In the next section we 
present a similar derivation for some cases of n = 3. 

7. n = 3 and cylindrical Bullough-Dodd 

The cylindrical Bullough-Dodd equation 

q"(t) + rV(t) = 4e 29 - 4e- 9 . (7.1) 

arises in the special case of n = 3 where c w = — u^c^-i. Then g x = 0, g 2 = — 93 and 
( [7. 1|) is satisfied by g = ^3. If we set £ := e 27rl//3 then may be chosen arbitrarily. If 
we choose it to be ((1 — £) then c^2 = C 2 (l — ( 2 ), c_i = and 

= sinvrs + 27r\/3 A sin(7r(s + 2)/3). 
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Again A is the one free parameter. The zeros given by 

1 3 • A x \ 
a ° = 4 ~~ 2tt arcsm 1 2 + 2\) ' ai ' " 2 = ~ a °' 

where A c = l/(2y / 3vr). Now A is the complement of the union of cuts 

(-oo, -3A C ] U [A c , oo). 

For A G A the zeros satisfy 

1 5 
Ra G (--, 1), »ai = l, »a 2 e(l, -), 

and so the extra condition is again automatically satisfied and our formulas hold for 
all A G A. If we write 

g(t)=Alog(^)-logfl + o(l) (7.2) 
then the connection formulas give in this case 

T f*±2\ T f2a±l\ 

A = -2a, B = 3- A ) 3 ( ) 3 ( , 

where we wrote a for a - F° r large t 

g(t) ~ 6Air 

Asymptotics at the critical value A c 

This section and the following ones are heuristic. Using the differential equation 
([7.1Q one can determine the correction terms to ( |7. 2| ) : 

, W = A ,og(I) + log B + ^ - ^ + ^ t -« 

(7.3) 

This is valid for A G A. To understand the higher order terms in more detail it is 
convenient to define 

w(t) = exp (-q(t)) 

where w satisfies the equation 

w " = -( W 'Y -l w ' + Aw 2 --. (7.4) 
w t w 
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The asymptotics we proved become the statement 

w(t) = Bt A (1 + o(l)) . 
Using (|7.4|) to calculate the higher order terms in the small t expansion for w we find 



24(A 2 -2A-2) , A 4 6 _ 3A 

+ (2 + A) 2 (l - A) 2 (A- A) 2 B + (1 - A) 4 (4 - A) 2 B 3 + " 



In contrast to (|7.3j ) the terms t 2m ~ 2mA only appear for m = 1 in the above expansion. 

As A varies from to A c a varies from to — |, so A varies from to 1 and B 
from 1 to oo. Observe that the first two terms in ( |7.5|) are of the same order in t as 
t — > (and A — > 1) whereas the others are of lower order. This suggests that when 
A = A c we have w(t) ~ tfli as t — ► 0, where 

/ t 2-2A \ 4 

fii := Jim 5 1- ^ 2(1 _ A)2 j = 21og(l/t) - - log 2 + 2 log 3 - 2j. (7.6) 



We now use the differential equation (7^) to find the higher order terms, which are 
polynomial in t and Q\. (The only property of f2i used in the formal expansion is 
dQi/dt = —2/t.) The expansion is 

W (t) = tn^-t 4 (nj + -fix + -)+^-t 7 (sin* + + 24 o^i + 8o)+o(t w n 4 1 ). 

9 V 3 9/ 2187 ^ ' 

(7.7) 

Thus (as for the n = 2 analogue [Q) if one were to alter the constant appearing in 
( [7.61 ) then the solution of ( |7.4| ) whose asymptotics is ( |7.7| ) would not match onto the 
solution that approaches 1 as t — > oo. 

These asymptotics at A = A c were checked by numerically solving ( |7.4p in both 
a foward and backward integration. There was agreement to nine decimal places at 
t = 1/4. 

Asymptotics at the critical value — 3A C 

We proceed as above and examine all terms that would be of the same order of 
magnitude as A — > — 3A C , when a — *■ 1 and A — > —2. These are the terms of the 
geometric series, those involving the powers + ^ A . Summing the series we see that 
we must compute 

Bt A 1 

hm 



*-i (1 - 2Bt 2+A / (2 + A) 2 ) 2t 2 (logt-log3 + 7 ) z 
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Defining 

Q 2 = log t - log 3 + 7 



we thus see that at A = — 3A 



c- 



1 



To compute higher order terms it is convenient to look at v(t) = l/w(t). Using the 
differential equation and only the property d^/dt = 1/t of ^2 we find 

v(t) = 2^ + ^ (sn- - |n£ + ^ - ^ + ^ - +o(t 14 ^ ). (7.9) 



Asymptotics for A > A c 

Think of A as being on the lower part of the cut [A c , oo). Then 

1 3 . 



a = — - 



where 



Thus A = 1 + 3ifi. Here again the first two terms in ( [7.51 ) are of the same order as 
t — > whereas the others are of lower order, and we obtain 



H := -arccoshfi + -^-), (/i > 0). 



w(t) = Bt A - 



t 



1-2A 



+ 0(t 4 ). 



B(l-A) 2 

Substituting in the values of A and B in terms of /i we find 



3 3i/i log(t/3) 



r(i/2-z>/2)r(-i/i) 

r(l/2 + i///2)r(i/x) 



If we had taken A to be on the upper part of the cut then we would have replaced 
/j, by —fi. The result would have been precisely the same. 

Remark. In [[| a method was described to find connection formulas for solutions 
of a class of equations including (|7.1|) . Away from the critical values the short- 
range asymptotics stated there correspond to the first two terms in (|7.5|) . As for the 
asymptotics at the critical values, our formulas agree with |J at A = — 3A C but at 
A = A c we differ by a factor of 2. 
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